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The semiclassical approximation to the coherent state propagator requires complex classical tra- 
jectories in order to satisfy the associated boundary conditions, but finding these trajectories in 
practice is a difficult task that may compromise the applicabihty of the approximation. In this 
work several approximations to the coherent state propagator are derived that make use only of real 
trajectories, which are easier to handle and have a more direct physical interpretation. It is verified 
in a particular example that these real trajectories approximations may have excellent accuracy. 

I. INTRODUCTION 

The path integral representation of the coherent state propagator K{zi, Z2,T) — {z2\e~^^'^/^\zi), 
where the \zi) are the usual harmonic oscillator coherent states, appeared in the works of Klauder 
and coUaboratorsii^ and of Weissman.^ A semiclassical, or stationary phase, approximation leads 
to classical trajectories satisfying Hamilton equations of motion but subject to special boundary 
conditions that can only be satisfied in a complex phase space. Aguiar and Baranger also considered 
this problem.4 g^^id discovered an extra term, which they called X, in the semiclassical approximation 
that had been overlooked in previous studies and that turns out to be essential for a correct theory- 
(the semiclassical spin propagator has a similar term^^ known as the Solari-Kochetov correction). 
Numerical calculations involving complex trajectories in the semiclassical coherent state propagator 
have been done for a variety of systems: Adachi considered a one-dimensional and time-dependent 
problem with chaotic dynamicsji Rubin and Klauder—, as well as Xavier and Aguiar— have treated 
ID bound systems; one dimensional tunnelling was considered in and also in 'lOl|; Van Voorhis 
and Heller presented calculations for one and two dimensioned^ and for the TV-dimensional Henon- 
Heiles potential; Ribeiro et al have worked with the 2D chaotic Nelson potentialji^ (numerical 
applications involving the spin coherent states have also appearedi^). Semiclassical approximations 
based on complex trajectories for the coordinate wave function, i.e. for the mixed representation 
K(x,z,T) = (a;|e~*^"^/''|z), were also developed, initially for the one-dimensional casei^^iSiii and 
then generalized to many dimensions«i^ The actual calculation of complex trajectories involves two 
difficulties: first, the effective dimensionality of the phase space is doubled, since both real and 
imaginary parts of position and momentum must be computed; second, the boundary conditions are 
defined part at initial time and part at final time, and finding the appropriate classical trajectory 
becomes a difficult problem known as 'root search'. Therefore approximations that make use only 
of real trajectories are certainly desirable. 

Since the propagator K{zi, z^, T) is a function of time, any complex trajectory that satisfies the 
boundary conditions at time T must belong to a whole 'branch' of trajectories, parametrized by T. 
In general, for a given system and for fixed values of zi and Z2 there are several such branches. In 
practice, once a solution is found for a particular value of T, one may obtain all elements of the 
same branch by making small steps forward or backward in time and using appropriate iteractive 
procedures. It may happen that for a certain value of time a relevant complex trajectory has a small 
(or even null) imaginary part, and in that case its branch was called 'nearly real' by Van Voorhis 
and Heller. dd'd^ It is possible that more than one 'nearly real' branch contribute to the semiclassical 
propagator for a given time, and thus one may consider only these branches and still accurately 
reproduce interference effects. 

A similar analysis can be made for the mixed propagator K{x, z, T), but in this case one usually 
holds T fixed and considers x as a parameter. Varying x thus produces a 'family' of trajectories, 
and again there may exist several such families. However, there is always a value of x for which 
the involved trajectory is real, and its family was called the 'main family' by Aguiar et al^ Using 
the main contribution alone is sometimes a very good approximation, but it can not reproduce 
interference because only one trajectory enters the calculation at a time. On the other hand, as 
already noted, finding all the necessary complex trajectories (i.e. performing the 'root search') is 
usually a difficult problem, specially in more than one dimension. Therefore the possibility was 
considerediSii^ of employing only real trajectories in the semiclassical approximation to K{x, z,T). 
This was done by approximating the complex trajectories by real ones, that are compatible with the 
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quantum uncertainties and satisfy less restrictive boundary conditions. The final real trajectories 
approximations are in principle less accurate than the original complex one, but they are much 
simpler and sometimes have practically the same accuracyiiSii^ 

The purpose of the present work is to present semiclassical approximations to K{zi, Z2,T) that 
are based only on real classical trajectories, thus making the calculation much more tractable. One 
method that accomplishes exactly this is the 'cellular dynamics', initially developed by Heller— 
(see also ^2^) and later generalized and applied to the stadium billiard with great successj^i This 
technique has shown to be accurate even for long times,^^'^'^ and it is actually very close in spirit to 
the present work, in the sense that the contribution of a complex classical trajectory is expanded to 
second order in the vicinity of a real one. However, Heller's starting point is the Van-Vleck-Gutzwiller 
formula for the semiclassical propagator, ^4 while we start from the formulation of Baranger et al^ 
and our results are slightly different from those of Heller. We also consider a variety of boundary 
conditions that the real trajectories may satisfy, something not discussed at length in,^^ 

Another approach to the semiclassical coherent state propagator that is based on real trajectories 
is the so-called Initial Value Representations, such as that of Herman and Kluk.^'* Recent reviews 
of this method can be found in |25j . Initial value methods are usually easy to apply and reasonably 
accurate for long times, but they require a numerical integration over all possible initial conditions. 
Since the present method requires only a few trajectories, at least for short times, it provides a much 
clearer physical picture. 

This article is divided as follows. In the next section we give a brief account of the semiclassical 
approximation to the coherent state propagator K{zi, Z2, T) and the complex trajectories. In section 
III we present the approximations that are based on the real trajectories defined by zi or by Z2. 
Real trajectories that satisfy mixed boundary conditions are investigated in section IV. We present 
an application to a nonlinear oscillator in sections V and VI and we conclude in section VII. 



II. THE SEMICLASSICAL COHERENT STATE PROPAGATOR 



The coherent states of a harmonic oscillator of mass m and angular frequency u> are defined by 

|z) = exp{za^ - z*a}|0), (2.1) 

where |0) is the oscillator ground state. The operators and a are respectively creation and 
annihilation operators, related to position Q and momentum P by 



The parameters b and c define natural scales of the problem, and are such that be — h and c/b — mu. 
It is easy to see that if we write 

then X and p are average values, 

{z\Q\z)^q, {z\P\z)^p. (2.4) 
The parameters b and c are related to quantum uncertainties, 

and we see that coherent states are minimum uncertainty states. 
These coherent states are never orthogonal, 

(Z2|zi) - exp - 1\Z2\^ + ^1^} , (2.6) 

and in the position representation they are Gaussians, 

{x\z) = n-h-^ exp j-^^^^ + ^p{x ~ q)| . (2.7) 
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In terms of the usual basis of number states |n), defined such that a'''a|n) = n|n), the coherent states 
may be written as 



-.7i. 



It is easy to see that they are eigenstates of the annihilation operator, a\z) = z\z). 

In order to write the semiclassical approximation to the quantum coherent state propagator 

K(zi,Z2,T) = (z2|e-*^^/'"'|zi), (2.9) 

we must consider a complex version of the phase space, i.e. we must make use of a coordinate q{t) 
and a momentum p{t) that are complex numbers. Following the approach of Ref. we define 

1 (q{t) p{t)\ 1 (q{t) .p{t)\ 

It is of fundamental importance to realize that v(t) is not the complex conjugate of u(t). In terms 
of these variables the boundary conditions become 

m(0) u' = zi, v{T)^v"^z*^. (2.11) 

There is nothing special about the values u{T) = u" and u(0) — u', they are to be determined 
dynamically. We use hereafter a prime (double prime) to denote initial (final) values, in order to 
simplify the formulas and stay close to the notation of 
The canonical coherent state propagator is 

i..(.„ T)^NY. exp { 1(5 + X)} , (2.12) 

where TV — exp{— ^jzip — ^|z2p} is a normalization factor, the summation is over all classical 
trajectories satisfying the boundary conditions, and the complex action is given by 



S{u',v'\T)^ / dt 
Jq 

This is related to the usual Hamilton action 



— (uu — vu) — Ti 



y(uV-uV'). (2.13) 



Sh^ f {pq~n)dt (2.14) 

JQ 



by 



S = Sh ^ y (w V - u V ). (2.15) 

The Hamiltonian that governs the classical movement according to the usual Hamilton equations 

■ 9n . an 

is the average value of the quantum Hamiltonian in coherent states, 

n^{z\H\z), (2.17) 

which is sometimes called the smoothed Hamiltonian. The quantity I is related to its second 
derivative, 

^-2 1 S**- <^'"*) 
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The prefactor in H2.12|l can be written only in terms of the complex tangent matrix. The clas- 
sical tangent matrix of a certain trajectory is the linear application that relates initial and final 
displacements about it. We take into account the quantum uncertainties to define it as follows: 



5q"/b\ _ ( \ fSq'/b 

Sp"/cJ \mpq TTippJ \5p'/c 

The complex tangent matrix, on the other hand, is defined as 

'Muu Muv\ ( Su' 



The relation between the matrix elements of these different representations is as follows: 



2M„„ 


= mqq + 


TTlpp 


+ inipq - 


ITflqp , 


2M™ 


= rriqq - 


TTlpp 


+ ITTlpq + 


ZTflqp , 




= rHqq - 


TTlpp 


ITTlpq 


ITflqp , 


2Myy 


= rriqq -+ 


TTlpp 


- inipq 4 


ITJlqp. 



It is possible to show that the second derivative of the complex action is given by 

i d^S 1 



h du'dv" ' 
and therefore the semiclassical coherent state propagator becomes 



(2.19) 



(2.20) 



(2.21) 
(2.22) 
(2.23) 
(2.24) 



(2.25) 



i^.c(.r,..T).g^exp{l(X + 5)^ (2.26) 



Upon fixing Zi and Z2-, the squared modulus of this propagator may be interpreted as a time 
dependent transition probability. On the other hand, if we fix zi and T and consider Z2 as a 
variable then \K{zi, Z2, r)P is a phase space representation, a Husimi function, of the evolved state 

If it happens that Myy tends to zero for a certain combination of [zi^ Z2,T), then we see that 
the semiclassical approximation (|2.26|) diverges. This is called a phase space caustio2iLi24i& and 
the quadratic approximation used in the derivation of (|2.26|) is not valid in its vicinity. In order 
to obtain an uniform approximation that remains valid at caustics it is necessary to employ a 
conjugate application of the Bargmann representation, as discussed in^^^^ We shall not be concerned 
with caustics in this work. 



III. THE 'LEAVING' AND THE 'ARRIVING' TRAJECTORIES 



We have seen that the classical trajectories entering the semiclassical propagator are determined by 
mixed boundary conditions. The initial position and momentum q' and p' are not the real numbers 
qx and pi, but rather some complex numbers such that u' = zi. Conversely, the final values 
are not 92,^*2 but are such that v" = It is in general not a easy task to find such trajectories 
in practice, even for simple systems. However, it may happen that the complex trajectory is close 
enough to a real one so that we may still obtain a reasonable result by expanding the propagator to 
second order in the vicinity of this real trajectoryiiSiiS We investigate this problem with some detail 
in the next sections. 



A. Leaving 

Let us suppose a certain complex classical trajectory that is to be used in the calculation of the 
semiclassical propagator, and let us assume it is not very different from the real trajectory that 
starts at the point (gi,pi). We call this the 'leaving' trajectory because it leaves the phase space 
point corresponding to the initial coherent state. After a time T the position and the momentum 
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will be some real numbers {qf,pf), generally different from the pair (172, ^2)- We will expand the 
complex action up to second order around this trajectory. If q' is the initial complex position and 
p' is the initial complex momentum, we may write 



q'^qi+Aqi, p'=pi + Api, 



(3.1) 



where Aqi and Api are assumed to be small (complex) quantities. Moreover, if q" is the final 
complex position and p" is the final complex momentum, we may write in a similar way 



q"=qf + Aqf, p"^pf + Apf. 
Therefore we have the approximation 



(3.2) 



S{u\v",T) W S{zi,Vr,T) + 



dS 



Aqi 



as 



Api 



1 

2 



Aqi 



dq'dp' 



AqiApi 



1 d^S 



2 dp'- 



ApI, 



(3.3) 



where the subscript r means that the quantity must be evaluated at the real trajectory (therefore 
v'^ = b~^qi — ic~^pi and u" = b^^Qf + '^^~^Pf)- order to obtain the derivatives of the action, we 
resort to equations p.lOfl and H2.15|l . Noticing that 



dS 



H 



OSh dq" dSn 



dq" dq' dp' 



OSh dq" 
dq" dp' ' 



one can obtain the derivatives of the total action, which are given by 



dS 



\/2 



dS_ 
d^ 



V2 



[v'^ - {rupp + imqp)u'l]. 



From the definition of the tangent matrix we have 



dq" 
97 



m 



99' 



dq" _ b dp" _ c dp" _ 



■ppi 



which determines, to first order, the final differences in terms of the initial ones: 

b 

-7 

c 

On the other hand, the boundary conditions 



Aq2 = mgqAqi + -viqpApi, Ap2 = -nipqAqi + rUppApi. 



q' , -P' qi , -Pi 
b c b c 



— - i— = — ~i — 
b c b c 



provide the secondary relations 

5-i[A(j2 + {qf - 92)] = ic-^[Ap2 + {pf -P2)], b-^Aqi = -ic-^Api. 
Solving for Aqi and Api in terms of (g/ — 52) and {pf — P2) we have 

Aqi = -Af^VK?/ - 52) - ib{pf - P2)/c], Api = icAqi/b. 
Substituting this in (|3.3|l one can see that the first order terms give 



dS 



dS 



g/ - 92 _ , Pf -P2 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



It is easy to take derivatives of equation H3.5|l in order to calculate the quadratic terms. In so 
doing we neglect derivatives of the tangent matrix elements, because this would be a higher order 
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correction. Adding up all quadratic terms and making the proper identifications, we see that it can 
be related to the difference (f" — zj) 

ih 

quadratic terms = -—MuvM~J{v'^ - z^ f. (3.12) 
Therefore the final result is the following: 

K,,P. {zi, Z2.T) = ^J=^ exp + Sr) + <K' - z^) + ^^K' - ^2)'} • (3.13) 

The subscript in Kq-^p-^ denotes that this formula was obtained using the 'leaving' trajectory. Notice 
that the prefactor and the extra term were not expanded but simply evaluated at the real trajectory, 
which is consistent with the original quadratic derivation of the semiclassical approximation. It is 
also important to remember that even though the action Sr is evaluated at a real trajectory, it 
continues to be a complex number. 

The expression ()3.13(l depends quadratically on the difference between the final value of the 
variable v along the real trajectory and the value that it would have in the complex trajectory. 
If by some reason the situation is such that and coincide, then this formula and the original 
one (|2.12|l will give the same result. One may argue that it is possible to obtain the same expression 
by expanding the action as 

{v';-zlf. (3.14) 

This is certainly true and actually an easy calculation. We have chosen the long way of using the 
position/momentum variables because this will be the only possibility in the next section. 



dS 



III * \ \ Q S 



B. Arriving 

What we call the 'arriving' trajectory is the real trajectory that starts in a certain initial point 
{qi,Pi) and after a time T arrives at the point (92, P2)- We can use this trajectory to approximate 
the semiclassical propagator in the very same way that we did with the 'leaving' trajectory. Similar 
to the previous arguments, we write 

q'^q, + Aqi, p' = p, + Ap^, q" = q^ + Aq2, p" = P2 + Ap2. (3.15) 

Inverting equation (|2.19|l we see that 

dq' dq' b dp' c dp' 



IP' a^// ~ h^P<i' a^ii ~ (3.16) 



Using these relations we can write the initial differences in terms of the final ones, analogously to 
what we did in (|3.7(l . Using the boundary conditions it is possible to show that 

Aq2 ^ -M-}[{q^ - qi) + i6(p, - pi)/c] , Ap2 = ~icAq2/b. (3.17) 

The first derivatives of the action arc in this case given by 

dS 



^ [u'; + [mpp - impg)v'^] , ^ 



- {mgq + imgp)v'^]. (3.18) 



We now expand the complex action to second order around this real trajectory. After simplifica- 
tions, we obtain 

Kg^P2{zi,Z2,T) = exp \ Ulr + Sr) + v'^{u'^ - zi) + l^P^iK - zi)A , (3.19) 

where the meaning of the subscript is evident. This time the expression depends on the difference 
between the initial value of the variable u in the real trajectory and the value that it would have in 
the complex one. Its interpretation is quite close to that of (|3.13|) . 
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IV. OTHER POSSIBLE REAL TRAJECTORIES 

In the previous section we saw that we may expand the semiclassical propagator in the vicinity 
of the real trajectories determined by the initial or by the final labels, {qi,pi) and (92,^2), which 
we called the 'leaving' and the 'arriving' trajectories respectively. Although these are probably 
the most natural real trajectories approximations, we can devise four more possibilities that are 
also interesting. Of course one may use any real trajectory to build an approximation -in fact, in 
principle it should be possible to find the 'best' choice by a variational approach, but this seems to 
be a highly nontrivial problem-, but the idea here is to obtain explicit formulas for the most natural 
cases. These are the four trajectories that are determined by pairwise combination of the coherent 
state labels. 

We shall present a detailed calculation for the case when the trajectories determined by the pair 
(91 5 92) are used. All other cases can be treated in a very similar way, and for them we shall be less 
explicit. 



A. From qi to 52 

Let us consider a trajectory which satisfies the following boundary conditions: it leaves qi at time 
zero and arrives at (72 at time T. Its initial and final momenta, pi and p/, remain unknown, but are 
real numbers. Differently from the previous section, now there may be more than one trajectory 
satisfying these requirements. We write 



q' = qi + Aqi, q" = q2 + Aq2, p' = p., + Api, p"=pf + Apf. 



(4.1) 



The initial and final momenta are regarded as functions of the initial and final positions. Therefore 
we may write 



di/_ 



Aqi 



dp' 



Aq2, Apf 



dp" 



Aqi 



dp" 



A92, 



(4.2) 



where again the subscript r means that the quantity must be evaluated at the real trajectory. On 
the other hand the boundary conditions u' = z\ and v" = imply that 



ic 



Api = -j^Aqi - {p. 



■Pi) 



Apf 



- — Aq2 ~ {pf -P2) 



(4.3) 



Since we are considering q' and q" as independent variables, the partial derivatives in (|4.2(l are 
given by 



dp^ 
dq' 



'99 



b m„ 



dp' 
'dq" 



c 1 

b TO, 



dp" 



c 1 
6 m„ 



dp" 
dq" 



(4.4) 



b m, 



where we have used that rnqqUipp 



1. Substituting this in 14.2|l and using (|4.3() we have 



Agi _ niqp [{pf - P2) - M2{pi - pi)] Aq2 _ rrigp [Mi{pf - P2) - {pi - Pi)] 
b ~ 



c 1 - M1M2 

where we have defined the complex numbers 

Ml = niqq + in 



1 - M1M2 



(4.5) 



(4.6) 



We now expand the complex action around this real trajectory up to second order, 

Aqi + -TT-fj Aq2 

dq" ^ 



1 9^ 

2 dq'^ 



dq'dq" 



AqiAq2 + - 



2 dq 



112 



Aqi 



(4.7) 
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Noticing that 



dq' 

we can obtain the derivatives of the total action, 



dSu 



Pf, (4.8) 



dS_ 



c dS 

-= [u'^ - Vlimqq + imqp)], — 

\ 2,m„r, Oq 



c 



V2r. 



[vl - u'J.{mpp + imqp)]. (4.9) 



After simpHfications, the linear terms can be written as 



linear terms = — ^[v'{pi — pi) — u"{pf — P2)]- (4-10) 
v2 

We now calculate the second derivatives and substitute H4.5II in H4.7f) . After many simplifications, 
the final result can be shown to be 

^9192(^1,^2, T) = Y] ff exp (^(X^ + Sr) + ^{Pf - P2) - ^{Pt -Pi)\ 
~: W(Myy)r I" v2c V2c J 



^ ''^p \ ~^(p' ~ pi)^ ~ T^^pf ~ ^^2)^ -^(p^- Pi)(pf - P2) , (4.11) 



where 



2 V 1 - M1M2 J 2 V 1 - M1M2 J 1 - M1M2 

This expression is more complicated that the ones we obtained in section III. This is so because 
the classical trajectories involved are determined by mixed boundary conditions, i.e. their initial 
and final positions. Its structure is nevertheless still the same: it depends on differences between 
the values of the variables in the real trajectories and the corresponding coherent state labels. The 
most important property of this formula is that the initial momentum pi is not known a priori. 
It must be determined as a function of the given parameters, and in fact there may be more than 
one possible value for it. Notice that since pi and also pf depend nontrivially on zi, Z2 and T this 
formula is not a simple Gaussian as it may seem at first. Once again, even though the function Sr 
is evaluated at a real classical trajectory, it will in general be a complex number. 

Notice that the differences pi — pi and pj — P2 are always divided by the momentum uncertainty 
c. Therefore only classical trajectories whose initial momentum is within a distance c from pi may 
be important for the semiclassical propagator. The same reasoning applies to the final momentum. 
We see that the real trajectories to be used in this formalism must be compatible with the quantum 
uncertainty principle. 

As a simple illustration of this formula, let us consider a harmonic oscillator of unit mass and 
angular frequency uj = c/b. An initial condition {q',p') leads, after a time T, to the final values 

q" = q' cos(wr) + — sin(wr), p" ^ ~ujq' sin{ijjT) + p' cos{ujT). (4.13) 

If we impose that the trajectory must start in qi and end in q2 then it is easy to see that there is 
only one possibility that satisfies these boundary conditions, for which 

_ U}{q2 - qi cosjuT)) _ uj{q2 cosjujT) ~ qi) 

sm((jjJ ) sm(cjJ ) 

In this case we have niqp = sin(a;r) and Mi — M2 — e'^'^^ , which leads to Ai ~ A2 ^ 1 and 
A12 = — e^*'^'^. The e*-^'/^ term cancels the prefactor. Finally, using that Sr = e^^"^"^ {qi/b + 
ipi/ c){q2/h — ipf /c)/2i we obtain 

if,,,,(zi,22,T) = exp{-i(|zi|2 + |z2|2) + e-'-rziZ2*}, (4.15) 

which is precisely the exact result. This comes as no surprise since the exact action in this case is of 
second order to begin with and thus all semiclassical approximations we consider in this work will 
be exact. 



9 



B. From qi to p2 

We now consider the real trajectory that starts in q' — qi with a certain momentum p' — pi and, 
after a time T, is in a final point q" = qf with the momentum p" = p2- We therefore treat q' and 
p" as independent variables, in which case we have the following partial derivatives: 



dp' c nipq dp' dq" 1 dq" b niqp 
dq' b TUpp ' dp" dq' nipp ' dp" c rupp 



(4.16) 



We may calculate the action's first derivatives, 

dS_ 
d^ 



d^ 

_ i^j, , ^j.y,^pp „,,i,pqyj, — — ■ 

V2mp„ dp" 



-[u" + v'^{mpp ~ impq)], 



[v'^-u'^{mpp + imgp)], (4.17) 



and after writing 

q' = qi + Aqi, q" = qf + Aqf, p' = p, + Ap„ p"^p2 + Ap2, (4.18) 
we may also obtain, using an expansion analogous to H4.2|) and the boundary conditions, the relations 

(4.19) 
(4.20) 



Aqi 


irupp [M2{pi - 


-Pi) 


- ic{qf - 


Q2)/b] 


b 


c 


1 + 


M2M3* 




Ap2 


irUpp [M^iqf 




- ib/c{pt 


-Pi)] 


c 


b 


1 + 


M2M* 





where M2 has already been defined and M3 = mpp + impq. 

After calculating the action's second derivatives, the final result is 

Kq,,A^uZ2,T) = E y=^-p |l(I. + Sr) - -|^(,,. - q2) - ^fe -Pi) 

X exp |-g(p, -pi)2 _ ^(^^ _ q,f + E^(p,^^p^)(q^ _ , (4.21) 

where the coefficients are given by 

^ _1 1 ( I-M2MA ^ _^ 1 f l~M;M* \ ^ impp 

^ 2 V 1 + M2M3* y ' ^ 2 V 1 + M2M3* y ' 1 + M2M3* ' ^ ' ' 

We see that the semiclassical propagator obtained is quite similar in structure to the one presented 
in the previous subsection. Only this time we have position and momentum in a more equal footing. 
As Pi — pi is always divided by c and qf — q2 is always divided by 6, we see that again the quantum 
uncertainties play a fundamental role in selecting the relevant classical trajectories. 



C. From pi to q2 



It is also possible to fix the initial momentum as pi and then search for an initial position q^ 
such that the final position is q2. In that case the final momentum will be some pf. Proceeding 
in complete analogy with the previous cases, we take p' and q" to be independent variables and 
calculate derivatives of q', p" and S with respect to them. After obtaining the values of Api and 
Aq2 in terms of {qi — qi) and (p/ — P2) and expanding the action to second order, the final result 
will be 

Kp,q^{zi,Z2,T) ^^^^=exp\Ulr + Sr) + ^{pf-p2)'' -^{q^ ^ qi) 

77 V(Mt,^)r I"- v2c V2b 

X exp|-^(g. - q,r - ^iPS-P^f - - 1i){Pf--P2)] , (4.23) 

where the coefficients are given by 

C, = l-ifi^^V C2 = l-iri^^V C,2 = ^^, (4.24) 
with A'h — m„„ + irririn- 
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D. From pi to p2 



Finally, we consider the trajectory determined by the pair {pi,p2). This has initial and final 
positions qi and g/, respectively. The procedure to obtain the semiclassical approximation is certainly 
clear by now, so it will not be repeated in any detail. The final result in this case will be 

where the coefficients are given by 



E. Summary of section IV 



In this section we have obtained four different semiclassical approximations to the quantum co- 
herent state propagator that are based only on real trajectories. The trajectories considered share 
the property that they are not determined by initial or final values, but satisfy mixed boundary 
conditions. Therefore finding them in practice is not trivial, but is certainly easier than finding the 
original complex ones. All the semiclassical propagators obtained are in principle able to reproduce 
quantum effects such as interference, since there may be more than one classical trajectory involved. 
They will be affected by caustics just like the original formula H2.26|l . but the location of such caustics 
will change because (M^t,)^. is different for each one of them. 

Which one of the several formulas obtained here and in section III is more accurate will depend on 
the particular problem at hand. We have considered only initial and final coherent states with the 
same value of the parameter 6, but a generalization of the semiclassical propagator was presented^ 
for more general 6's, and the calculations presented here may be adapted to that case with no 
essential difficulty. Let us suppose for a moment that the initial coherent state \zi) has a position 
uncertainty bi while \z2) has a position uncertainty 62- If these numbers are small that means the 
states are very narrow in the position representation, while having a large uncertainty in momentum. 
In that case we conjecture that an approximation in the spirit of Kq-^q^ would be the most effective 
one. If bi is small but 62 is large, than Kq-^p^ would be a better candidate, and so on. Of course for 
the free particle and the harmonic oscillator they are all exact, regardless of the values of 61 and 62. 

In the next section, we present an application of the formalism just presented to a nonlinear 
system. The purpose is not to attempt an exhaustive investigation of the several possibilities, but 
rather to illustrate the method in a simple case. We have chosen a system for which many analytical 
results are possible so that the main properties of the theory do not disappear under numerical 
calculations. 



V. APPLICATION TO A NONLINEAR OSCILLATOR: SHORT TIME 

We consider the nonlinear Hamiltonian 

H^n.ia^ar^pJ^±^^^^^^^, (5.1) 

which is diagonal in the usual number basis, 

H\n) = En\n) = hujn^\n). (5.2) 
The quantum propagator in this case is quite simple: 

K{z,,Z2,T) = {z2\e-^"^/^\z,) = V iflflLe--^-^. (5.3) 

■^-^ nl 

n=0 
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We shall be interested, for simplicity, only in the diagonal case 

oc I |2n 

K{z,,z^,T) = e-l-^l^ ^e-»^-^, (5.4) 
^ — ' n! 

n=0 

whose squared modulus is the return probability, 

p(zi,r) = |x(zi,zi,r)|2. (5.5) 

This function is periodic with period Tr = 2tt/lu. In the semiclassical limit the term that is respon- 
sible for the largest contribution to the sum in 1)5. 4|l is no « jzip. If we linearize the exponent in 
the vicinity of this term we have 

2 



P(zi,T)«e-2|^i 



(5.6) 



Notice that expression H5.6|l has a distinct time scale. 



Tc^—. (5.7) 
noui 

The quantities Tr and Tc are usually called revival time and classical timCfSS, 

For short times, we can approximate \K{zi,T)\'^ w 1 — {{H^) — {H)^) T^/fi^, where (•) denotes 
an average value in the state \zi). For the system in question, this gives 

F(zi, T) « 1 - (4|zi|6 + 6|zi|4 + |zi|2) cj^T\ (5.8) 

Let us write zi — {q + ip)/^/2 and take for simplicity the value q ^ 0. Since the movement in phase 
space has circular symmetry, this choice is of no fundamental importance. The short-time expansion 
1)5. 8|l becomes simply 

P(;2i,T)«l-i(/ + 3/+p2)T2. (5.9) 

Let us now turn to the semiclassical approximation. From now on we set h — lu ^ 1, which implies 
b = c = 1. The smoothed Hamiltonian associated with 15. 1|) is 

(p^ + g^)(p^ + g^ + 2) 
n — =uv[uv + l), (5.10) 

and the corresponding Hamilton equations are 

. dn . dn 

q= — = ap, p = = -(T(j, (5.11) 

where we have defined a = p^ + (f' + 1. If we note that {cr, 7i} = 0, and thus that cr is a constant of 
the motion, then it is clear that 

q" = q' cos{ijt) + p' s\\\{at), p" ^ p' cos{at) - q' sm{crt) . (5-12) 

We see that the classical trajectories have a period of motion that is energy-dependent and given by 
27r/(T. If we remember that Uq -I- 1/2 = (5^ -)- p^)/2 we see that in the semiclassical limit this time 
scale becomes precisely Tc- 

The tangent matrix that is associated with the classical trajectory that starts in {q',p') can be 
obtained by simply differentiating the equations of motion. We must remember that the angular 
frequency a is not uniform. The result is 



cos(cri) sm{at)\ fl + 2q'p't 2p'^t 



^iripq nippj sin(o-t) cos{at) J —2q'^t \ — 2q'p'tJ ' (5.13) 

The action of such a trajectory is easily seen to be 

{a-lfT . (g-l) 
S= « — ^ — , (5.14) 

while the extra term is 

X={a-l/2)T. (5.15) 

The result of this semiclassical approximation based on complex trajectories will be given in section 
VI. 
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A. The 'leaving' and the 'arriving' trajectories 

The first possibility we consider is to approximate the return probability by using only the real 
trajectory that leaves the position q = with momentum p. The tangent matrix for that trajectory 
is 

_ f cos(at) sm{c7t)\ (I 2p^T\ 
siniat) cos{at) J \0 1 J ' 

which gives the values 

= -ip^Te'"'^, - (1 + ip''T)e"'^. (5.17) 

The angular frequency is ti = + 1 and the final points in phase space are g/ = psin(crT) and 
Pf = pcos(crr). This corresponds to 

< = {v':r = i^^-'^^- (5-18) 

Inserting all this information, together with H5.14|l and (|5.15|) . into the formula H3.13|l we have 
Kq,p, (zi, T) = \ exp ( f (/ + 2p') - sin(ar/2) 

The first thing we note is that for p = we obtain the exact result K = 1. Moreover, in the short 
time regime we can expand (15.19(1 and obtain 



|X,,p,(2i,T)|' w 1- - (p6 + 3p4+p2^r2 (short times), (5.20) 



which again reproduces the exact calculation. For later times we should not expect exact agreement. 

Let us now turn to the 'arriving' trajectory, the one that starts in qi,Pi and arrives at te position 
g = with momentum p after a time T. The tangent matrix is less trivial than in the previous case, 
but in the end we get 

M„„ = -ip'^Te-'"'^ , M^y = (1 + ip^T)e"'^. (5.21) 

Since + is a conserved quantity, we have a = qj^ + pj + 1 — p^ + \ . After the whole calculation 
is done, we find out that Kq^p^ — Kq-^p^ . This indicates that perhaps these two approximations will 
always have the same content of information, something that is not completely unexpected because 
of the dual role of \zi) and \z2) in the quantum propagator. 



B. The qi q2 possibility 

In that case the trajectories that enter the approximation have initial momentum given by the 
equation 

p, sin[(p2 + i)T] ^ 0. (5.22) 
Of course one solution to this equation is 

p,{=Pf)^0, (5.23) 

in which case the particle simply stay still and a = 1. It is easy to see that for this trajectory the 
tangent matrix is very simple, 

cos(T) sin(T)\ 
-sin(r) cos(T) j ' 



fmqq 


mqp\ 


-{ 


\mpq 


Vflpp J 
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and therefore Mi = = M^^ = e*-^. The contribution of this trajectory to the propagator is 

Kq = exp{-ip2 sin(r/2)e"*^/2}, (5.25) 

where we have used Sr+Ir — T/2. Notice that for p = we again have the exact result — I- We 
also note that the function \Kq\'^ has a period of 2-k, which of course corresponds to the quantum 
revival time. 

We now turn to the other solutions of equation (|5.22|) . They are of the form 



pUn) = — - 1' (5-26) 

which leads to — 2n7r/T. In this case we have less trivial trajectories, for which the tangent 
matrix is given by 



Ipq Tflpp 



1 2piT 
1 



(5.27) 



and we see that the prefactor is = (1 + ipfT), while Mi = M2 = —ApfT'^. The action and the 
extra term are given by Sr = {pfT — 2ip|)/4 and Tr — {pf + 1/2)T. The coefficients in H4.11|l are 

^1 = ^2 = 1 + „./5^2T-V ^12 = -—^. (5.28) 
2(1 + ipfT) 1 + ipfT 



After many simplifications, we obtain 



Tl = 



where the contribution of the trajectory with label n (different from zero) is given by 

if„ = \ exp{--^4^(K - pf + '^{pt + ipl + 2)}. (5.30) 
y/l + ipfT 1 + ipfT 4 

Note that for short times Pi{n) is very large, so Kn becomes negligible and Kq gives the only 
contribution. However, it predicts the initial decay \Kq\'^ w 1 — 2p^T^, which is very slow compared 
to the exact calculation H5.9|l . The two results agree only for very small values of the momentum p. 

Concerning the contributions Kn-, we see that for a given instant of time the value of n that will 
contribute the most is that for which Pi{n) is as close as possible to p, because of the Gaussian 
decay in (|5.3U|) . If we impose pf{n) « p'^ we have T « 2n7r/(p^ + 1), which means that the return 
probability has a maximum at the classical period, in agreement with the exact result. 



C. The qi p2 possibility 

If we impose that the classical trajectory starts in g = with momentum pi and ends at with 
momentum p, we have 

g/ = Pi sin (o-T), p = piCOs{aT), a=p1 + l. (5.31) 

These transcendental equations have no explicit solution. If we confine ourselves to the short time 
regime, then we can write pi ^ p and qf w p{p^ + 1)T. The tangent matrix is given by 

rugq ■mqp\ ^ f cos(a-t) sin(a-t)\ fl 2pfT\ (5 32) 

nipq ruppj \-- sin(CTt) cos(crt) y 1^0 1 y ' \ ■ ) 

and we obtain M2 = e*'^^(l + 2ipjT) and A/3 = e~*'^^ - 2plT. Substituting this in (|OT|l . we obtain 
\Kq,p^'^ « 1 - (p^ + V +p2 J y2 ^gj^Qj.^ iijnes), (5.33) 



which decays faster than the exact result but is a better approximation than the one obtained using 
the qi — > q2 trajectory. We see that the different approximations may lead to very different results. 
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D. The pi q2 possibility 

The equations of motion in this case are 

= q,cos{aT) +psm{aT), pf ^ ~qiSm{aT) + pcos{aT), a = q'^ + p'^ + 1, (5.34) 
while the tangent matrix is 



^qq '"'qp 
Ipq T^pp 



cos(cri) sm{(7t)\ C 1 + 2qipT 2p^T 



(5.35) 



The situation here regarding solubiUty of the equations is even worse than in the previous case. 
Once again we restrict the analysis to the short time regime. Then it is possible to write the first 
equation as qi ~ —paT/2 and find 



1 - v/l-4p2T2(p2 + i)) ^ (5.36) 



* 2pT 

which we substitute in the first equation to find pf. Carrying out the whole calculation will give in 
the end 

\Kp,q2 P ~ 1 - ^ + Sp" + P^) (short times), (5.37) 
which agrees with the exact result. 



E. The pi p2 possibility 

Finally, in the last possibility we have 

g/ = cos((tT) +psin(crT), p = — sin(crr) +pcos((7T), a — qf+p'^ + l. (5.38) 
In the short time limit we have again 

9* « (l - Vl- Vr2(p2 + 1)) , (5.39) 

and the final result is 

\Kp^P2? ~ 1 - y?"^ (short times). (5.40) 

This is kind of intermediate between the result we found in subsection B and that of subsections C 
and D. 

VI. APPLICATION TO A NONLINEAR OSCILLATOR: NUMERICAL RESULTS 

Before we consider the semiclassical approximations based on real trajectories for longer times, 
let us see how well the original one (|2.26|) compares to the exact result. This has been considered 
in detail in poj . so we just present the result. Given the initial condition u' = zi, for each time 
T we must find a value for v' such that v" = Z2 = z^. This problem usually has more than one 
solution, and we must add their contributions coherently. In Fig.l we see the return probability 
as a function of time (in units of Tc) for the case p = 10, which we have chosen to ensure that we 
are in the semiclassical limit. The corresponding classical period is Tc ~ 0.062. The exact and the 
semiclassical results are indistinguishable in this scale. 

In the previous section we saw how the different approximations based on real trajectories per- 
formed in the short time regime. The exact result was reproduced only by the 'leaving' and the 
'arriving' formulas and by Kp-^q^. We now turn to the less simple case of arbitrary T, when the 
classical trajectories and the associated propagators must be computed numerically. 
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Let us start with the propagator Kq^p^ , which is based on a real periodic orbit. Its initial decay is 
exact, and we can see how well it does for later times in Fig. 2. It is able to reproduce the height of the 
peaks with great accuracy, but not their widths. Since there are never more than one contribution 
for each time, it never displays any interference effects. 

This is not the case for Kq^q^. We see from (|5.29|) and (|5.30l) that it consists in the sum of many 
contributions. We focus on the values n = 1,2,3. Their individual contributions are depicted in 
Fig. 3. Notice that the second and third peaks overlap. When we calculate the total propagator, 
this gives rise to interference. The final result is indistinguishable from the exact one (for T > Tc/2, 
because we have ignored Kq which gives a bad initial decay). 

So far the propagators could be obtained analytically. Since the calculation of Kq^p^ depends on 
the solution of the transcendental equation (|5.31l) . we must resort to numerical routines. Let us try 
to find solutions to the second equation in H5.31|l in the vicinity of the first period, T « Tc- In Fig. 4a 
we see that there are two solutions (solid lines) for T <Tc and no solution at all for T > Tc- This is 
because the cosine function with a real argument is always less than unity, and thus Pi must always 
be greater than p. The complex solutions do not have this obstruction, as we can also see in Fig. 4a 
(dashed line), where we plot the real part of the complex momentum that satisfies the boundary 
conditions H2.11|l . Therefore the semiclassical approximation based on the complex trajectory can 
reproduce the whole peak, while Kq-^p^ is discontinuous. 

In Fig. 4b we see the squared modulus of the exact propagator and the values of {Kq-^p^ P obtained 
using the two available real trajectories. Note that one should not add these results. They are 
independent and we may choose any of them, because both real trajectories are good approximations 
to the actua l com plex one (the real trajectories do not come from a saddle point approximation). As 
observed in [Ml, the mixed propagator (x|e *^'^/''|z) can also be discontinuous when calculated 
using real trajectories. But in that case there are caustics involved, while here we have an algebraic 
obstruction. 

The discussion of the approximation Kp-^q^ is quite similar to the above. The solutions to the first 
equation in (|5.34|) arc shown in Fig. 5a, where we also show the real part of the complex position that 
satisfies the boundary conditions (|2.11() . Again there is no solution for T > Tc and the semiclassical 
propagator is discontinuous, as we appreciate in Fig. 5b. The results are practically the same as in 
Fig. 4b. 

Finally, the propagator Kp-^p^. This time we solve numerically the conditions H5.38|l and find that 
there is a single real trajectory for T < Tc and no one for T > Tc- The final result is in Fig. 6. 

VII. CONCLUSIONS 

Several approximations to the semiclassical coherent state propagator {z2\e~^'^'^^^\zi) were pre- 
sented that are based solely on real classical trajectories. Two of these approximations do not involve 
mixed boundary conditions and thus are not hindered by the associated 'root search' problem. The 
remaining four possibilities are based on trajectories that are determined by initial and final data, 
but since they are real for all times they are simpler to determine than the original complex ones. 

As a testing ground we have used the nonlinear system H = (a^a)^. Only one of the approxi- 
mations, namely Kq-^q^, reproduced the exact result to the fine details. This is certainly due to the 
particular initial coherent state that was chosen, one corresponding to g = and p — 10. Had we 
chosen for example g = 10 and p = and then Kp^p^ would give excellent results. We could also 
consider a nondiagonal propagator, and in that case we would expect Kq-^p^ , for example, to improve 
its performance. 

Straightforward extensions of this work include the already mentioned case of different position 
uncertainties (squeezed states) and also higher dimensional systems. It is also possible to fix the 
time T and the initial state \zi) and to regard \K{zi, Z2,T)\'^ as a Husimi function defined in the 
Z2 plane. This is technically more difficult than what we have presented here, because it involves 
finding classical trajectories -usually more than one- parametrized by points in the plane. 

Similar results can be obtained for the semiclassical SU{2), or spin, coherent state propagator, 
even though the introduction of position and momentum variables in that case is not as natural. 
The associated phase space is also two-dimensional, but since it has curvature the calculations may 
be a little more involved. The same may be said about the semiclassical SU{1, 1) coherent state 
propagator. Since these groups have wide applications, it would be interesting to also have the 
corresponding approximations based on real trajectories. 
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FIG. 1: Squared modulus of the exact propagator for g = and p = 10. The semiclassical approximation 
based on complex trajectories is indistinguishable from it in this scale. Time is in units of the classical 
period. 




FIG. 2: The function ^ as a function of time. It reproduces well the height of the peaks, but not their 

widths, and shows no interference. 




FIG. 3: The individual contributions \Kn\^ for the approximation Kq^q^- We show the cases n = 1, 2 and 3. 
When they are added there is interference, and the exact result of Fig.l is reproduced with extraordinary 
accuracy for T > Tc/2 (we have not included Kq in the calculation). 
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T 

FIG. 4: Top: real solutions to the equation p = pi cos((pf + 1)T) in the vicinity of the classical period (dashed 
lines). Wc also show the real part of the momentum for the complex trajectory (solid line). Bottom: 
approximation IJfgipsI^ (dashed lines) compared to the exact result (solid line). Since there are no real 
trajectories for T > Tc, the propagator becomes truncated. For T <Tc there are two possibilities. 
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FIG. 5: Top: real solutions to the equation = qi cos(aT) +pcos{aT), where a = g| + 1, in the vicinity 
of the classical period (dashed linos). Wo also show the real part of the position for the complex trajectory 
(solid line). Bottom: approximation I^J^piqal^ (dashed lines) compared to the exact result (solid line). The 
situation is analogous to the previous figure. 
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FIG. 6: Approximation ji^pipsl^ (dashed line) compared to the exact result (solid line). This time only one 
real trajectory exists for T < Tc, but again the propagator is truncated. 



